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ABSTRACT 



Aims. We present the first full-star three-dimensional explosion simulations of thermonuclear supernovae including parameterized 
deflagration-to-detonation transitions that occur once the flame enters the distributed burning regime. 

Methods. Treating the propagation of both the deflagration and the detonation waves in a common front-tracking approach, the 
detonation is prevented from crossing ash regions. 

Results. Our criterion triggers the detonation wave at the outer edge of the deflagration flame and consequently it has to sweep 
around the complex structure and to compete with expansion. Despite the impeded detonation propagation, the obtained explosions 
show reasonable agreement with global quantities of observed Type la supernovae. By igniting the flame in different numbers of 
kernels around the center of the exploding white dwarf, we set up three different models shifting the emphasis from the deflagration 
phase to the detonation phase. The resulting explosion energies and iron group element productions cover a large part of the diversity 
of Type la supernovae. 

Conclusions. Flame-driven deflagration-to-detonation transitions, if hypothetical, remain a possibility deserving further investigation. 
Key words. Stars: supernovae: general - Hydrodynamics - Instabilities - Turbulence - Methods: numerical 



1. Introduction 

The observables of Type la supernovae (SNe la) are determined 
by the composition of the ejecta and the energy release in the 
explosion. Both depend sensitively on the density at which the 
carbon-oxygen white dwarf (WD) material is burned by the ther- 
monuclear flame. This in turn is a consequence of the competi- 
tion between the flame propagation and the expansion of the WD 
due to nuclear energy release. Therefore, reproducing the char- 
acteristics of SNe la in a thermonuclear white dwarf explosion 
model requires a specific sequence of flame velocities. 

After the thermonuclear flame has ignited near the center 
of the WD, two distinct modes of outward propagation are hy- 
drodynamically admissible. In a subsonic deflagration the flame 
is mediated by the thermal conduction of the degenerate elec- 
trons, while a supersonic detonation is driven by shock waves . 
Obser vational constraints rule out a prompt detonation (lArnetlL 
Il969t> and thus the flame has to start out as a deflagration 
(Nomoto et al.LI 1976h . Such a scenario leads to a stratification of 
dense fuel atop light ashes. Due to Rayleigh-Taylor and Kelvin- 
Helmholtz instabilities, turbulence is generated. This results in 
a turbulent cascade with eddies decaying into ever smaller ones 
until they are dissipated at the Kolmogorov scale (< 1 mm). In 
the cascade, the turbulent velocity decreases with length scale. 
Thus, turbulent motions interact with the flame down to the so- 
called Gibson length lc, where the turbulent velocity fluctuations 
v' become comparable to the laminar flame speed s\, 

v'Qg) ~ si- 



For most parts of the explosion, lc is much larger than the 
flame width If, indicating that the internal flame structure is un- 
affected by turbulence. Hence the interaction in this flamelet 
regime of turbulent combustion is purely kinematic, corrugat- 
ing the flame and accelerating it by increasing its surface. 
However, as the WD star expands, the fuel density dr ops, s\ de- 
creases, and Zf increases (Timmes & Wooslev, 1992). As soon 
as Iq <, If, turbulent eddies distort the internal flame struc- 
ture. This distributed burning regime commences when the fuel 
density has dropped close to values at which deflagration burn- 
ing ceases dNi emever & Wooslev, 19971: [Niemever & Kersteinl 
1997l lRopke & Hillebrandtll2005bh . 
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In order to burn significant fractions of the WD, 
a high flame propagation speed is necessary. Even a 
turbulence-boosted flame is relatively slow so that it may 
not lead to a sufficient fuel consumption before fuel 
densities drop below the burning th reshold (<10 7 gcm~ 3 ) . 
Although turbulent deflagration models dReinecke et all 120 02b ; 
iGamezo et all 120031: iRooke & Hillebrandll l2005at iRopke et al., 
2006; ISchmidt & Niemeyerl 120061) reproduce the observ ational 
features of weak SN la explosions (Blinnik ov et al.L 12006). there 
exist open questions dKozma et all 20051) and presently it seems 
that they cannot account for the more energetic events. This issue 
certainly needs more exploration (in particular the buoyancy in- 
duced turbulence effects at the high Reynolds n umbers of ~10 14 
expected in SNe la, see ICabot & CookL 120061) but it may also 
be interpreted as an incompleteness of such models. If more 
burning is necessary than the deflagration model allows for, a 
deflagration-to-detonation transition (DDT) would provide the 
ultimate flame a cceleration . The problem of such delayed d et- 
onation models (Khokhiovl Il99lt IWoosley & Weaverl 11994 is 
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the absence of a robust physical mechanism that allo ws predict- 
ing a DDT in thermonucl ear supernova explosions (Nie meveri 
1999: iRopke et all l2004allbh - Therefore invoking a DDT makes 
it necessary to introduce a free parameter into the model. 

Exa mples of suc h parametrized models were recently 
given. Game zo et al. I (12005b assume the DDT to take 
place at an arbitrarily chos en instant and location, while 
iGolombek & Niemeveri (120051) impose a DDT once the flame 
enters the distributed burning regime. 

Here we present the first three dimensional simulations 
of delaye d detonations based on t he num erical methods pro- 
posed by Golombek & Niemever (120051) and adopting their 
DDT parametrization. Since the only instant where the flame 
properties change drastically is the beginning of the distributed 
burning regime, it seems natural to assume a DDT here 
dNiem ever & Woosl eylll997l) . Although presently there is no in- 
dication of a DDT occurring in the distributed burning regime 
dLisewski et all l2000t IZingale et all [2005). more exploration is 
necessary to settle this question. Our DDT parametrization does 
not a priori fix the instant of its occurrence nor its location on the 
flame and thus a failure triggering a viable detonation is possi- 
ble. If successful, however, the detonation propagation is likely 
to proceed asymmetrically. In order to avoid artificial symme- 
tries in the setup of the simulations, we pursue a full-star imple- 
mentation of the model. Our aim is to determine whether such 
a model is consistent with global properties of observed SNe la, 
which is non-trivial since here the DDT is not used to tune the 
model but results from a physically motivated hypothesis. 

2. Modeling approach 

Aiming at parameterizing the location and time of the DDT ac- 
cording to properties of the deflagration flame, a self-consistent 
treatment of the interaction of the flame with turbulence is essen- 
tial. This is achieved by a large -eddy simulation approach on th e 
basis of metho d s pres e nted by Niemeyer & Hillebrandtl (|1995 | ) . 
IReinecke et ail d 19991), IReinecke et alj d2002al) . iRopkel {2005), 
and lSchmidt et alj d2006l) . In particular, the deflagration flame is 
represented as a sharp interface between fuel and ashes applying 
the level set method and unresolved turbulence is described via a 
subgrid-scale model from which the turbulent flame propagation 
spe ed is determined. 

iMaier & Niemeveri ((2006) found that even tiny regions of 
ash stall the detonation wave. Consequently, the detonation has 
to sweep around the complex deflagration structure. This re- 
quires a versatile i mplementation of the detonat ion. Following 
the suggestion of Golombek & N iemeveri d2005l) we model the 
detonation wave in the same approach as the deflagration flame 
applying a second level set treatment (see this reference for thor- 
ough tests of the implementation). 

Lacking a sound theoretical understanding of the mecha- 
nism of the DDT, we restrict our model to a single DDT which 
is achieved by initiating the detonation level set the first patch 
of the deflagration flame where If > 0.98 Zg- Although this ap- 
proach may seem rather artificial, it was chosen to simplify mat- 
ters and to obtain a clear picture of the detonation propagation 
through the WD. Physically, it is motivated by the assumption 
that a DDT is a rare event. 

In fuel regions the propa gation spee d of the detonation front 
is prescribed according to Sharpe ( 1999), while propagation into 
ash is suppressed. One question to be addressed in the present 
study is to which degree such a handicapped detonation is capa- 
ble of burning major parts of the WD matter - an issue calling 
for three-dimensional simulations. 



Our description o f the nuclear reactions is based on 
IReinecke et al.l d2002ah with the exception that the lo wer thresh- 
olds ar e set differently for burning C and O, following lKhokhlovl 
(1991). While we allow C consumption down to fuel densi- 
ties of 10 6 g cirT 3 , O consumption ceases already at 10 7 g cm -3 . 
The WD is constructed using a suitable equation of state as- 
suming isothermality at T - 5 x 10 5 K, a central density of 
2.9 x 10 9 gcirT 3 , and a composition of equal parts by mass of 
C and O. It is set up on a Cartesian computational grid of 256 3 
cel ls with outflow boun dary conditions on all sides. As described 
by IRopke et alj d2006l) . we follow the propagation of the defla- 
gration flame and the WD expansion with two nested moving 
grid patches optimizing the flame resolution. 



3. Simulations 

With the DDT parametrization fixed as above, the only major 
undetermined parameter that impacts the evolution of the model 
is the configuration of the deflagration flame ignition. Although 
not attempting a thorough exploration of the parameter space 
here, we set up three different models in the multi-spot sc enario 
similar to the simulations presented by Rop ke et ail (|2006). This 
allows us to test the robustness of the imposed DDT criterion. 
By placing 5 (Model DD_005), 20 (Model DD.020), and 800 
(Model DD_800) flame ignition kernels inside a radius of 1 80 km 
around the WD's center, the efficiency of burning in the deflagra- 
tion phase is significantly altered. While Model DDJ305 leads 
to relatively little burning, igniting 800 flames results in a strong 
deflagration phase. We focus our description on the intermediate 
case of Model DD_020. 

After ignition, the deflagration phase p roceeds in a simi- 
lar way as described by Rop ke et all |2006). The burning bub- 
bles grow by flame propagation and rise buoyantly towards 
the surface of the WD. Due to instabilities and partial merger 
of the bubbles, a complex connected structure develops. The 
parametrized DDT criterion is met first 0.724 s after ignition at 
the outer edge of the flame front. Here, the detonation is trig- 
gered by initiating the corresponding level set, as shown in the 
upper panel of Fig. [T] The fact that the outer edge of the de- 
flagration flame is chosen by our DDT criterion is not surpris- 
ing. Turbulence is generated preferentially at large buoyant bub- 
bles and at the same time the density is lowest at the outer- 
most parts of the flame making the flame broadest here. This 
favors the transition to the distributed burning regime and thus 
the parametrized initiation of the detonation wave. 

The center and lower panels of Fig. [1] show the propagation 
of the detonation wave. As implemented in our model, it can- 
not cross ash regions and, burning towards the star's center, it 
therefore wraps around the corrugated deflagration structure. In 
this way, it takes about 0.2 s before it arrives at the center of 
the WD. Meanwhile, the star keeps expanding and the deflagra- 
tion continues in regions not yet reached by the detonation wave. 
Consequently, the density of the fuel ahead of the detonation 
drops quickly after passing the WD's center. In Model DD_020 
burning stalls shortly before the detonation reaches the far side 
of the deflagration structure. This still implies burning of most of 
the WD, since these deflagration features have already reached 
the low density edge of the star. 

The DDT parameters and the results of all three simulations 
simulations are summarized in Table Q] The more kernels the 
deflagration is ignited in, the less energetic becomes the explo- 
sion and the less NSE elements are synthesized. The latter effect 
is partially compensated by more IME produced in weaker ex- 
plosions giving rise to comparable masses of total burned mate- 
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Fig. 1. Initiation and propagation of the deflagration wave (white 
isosurface) in Model DD_020. The deflagration flame is shown 
as dark isosurface and the extend of the star is indicated by the 
central plane mapping the logarithm of the density. 



rial. In terms of the masses of the explosion products, the entire 
fuel remaining from the deflagration is affected by the detona- 
tion in all simulations. This is the reason for the virtually com- 
plete C-depletion in all models. The different expansion history, 
however, is reflected by up to ~ 0.2 M Q of O originating from 
material detonating at densities less than 10 7 g cirT 3 . 

Although the fuel density at the location of the DDT in- 
creases with the number of deflagration ignition spots and the 



DDT sets in earlier, the explosions do not become stronger. The 
most important parameter for the overall explosion strength is 
the nuclear energy release in the deflagration phase prior to the 
DDT, £™ T . It accounts for the expansion of the star and there- 
fore determines which fraction of fuel is available for the deto- 
nation to burn to NSE or IME. Obviously, the asymptotic kinetic 
energy is larger in models releasing less nuclear energy prior the 
the DDT (cf. Table [TJ. Moreover, the complexity of the defla- 
gration flame increases with more ignition spots impeding and 
delaying the propagation of the detonation. These effects cause 
a larger number of ignition spots to shift the emphasis of the 
explosion model from the detonation phase to the deflagration 
phase which is characterized by the ratio of the mass of nuclear 
ash produced by the deflagration to that produced by the detona- 
tion, ftdtf (see Table[U- Whereas in Model DD_800 the deflagra- 
tion mass consumption exceeds that of the detonation by a factor 
of 1 .59, the detonation burns about three times as much material 
as the deflagration in Model DD_005. 

This way, we obtain three different final configurations. 
Assuming that a large part of the NSE material is actually 56 Ni, 
which by radioactive decay powers the light curve of the event, 
Model DD_005 would correspond to a very bright supernova. 
Model DD_020 provides an intermediate case with values sim - 
ilar to the one-dimensional W7 model ( Nom oto et all [l984), 
which is considered to be prototypical for many SNe la, while 
Model DD_800 corresponds to a dimmer event. In the latter, the 
detonation does not significantly enha nce the NSE produ ction 
over similar pure deflagration models (Ropke et al., 2006). The 
major difference to pure deflagration models is a more complete 
burning of the outer parts of the WD leading to larger asymp- 
totic kinetic energies of the ejecta and a different composition 
of these. The inner core is dominated by NSE material, par- 
tially mixed with IME originating from the deflagration phase 
and low-density detonation burning. Above, a layer dominated 
by IME is found which is completely C-d epleted in agreeme nt 
with a recent analysis of observed spectra (Steh le et al. , l2005l) . 

The asymptotic kinetic energies of the ejecta range around 
the expected value of ~10 51 erg. In all tested scenarios, the DDT 
occurs before the deflagration structure encloses large pockets 
of fuel, which would be unreachable for the detonation wave. 
This, however, may in part be due to the low resolution of the 
simulations and needs further exploration. 

4. Conclusions 

In the presented study, we imposed a hypothetical DDT criterion 
on SN la explosion simulations. Being guided by a physical sce- 
nario (the transition of the deflagration flame from the flamelet 
to the distributed burning regime), it cannot be used to tune the 
models. Thus, the success of the presented simulations to repro- 
duce gross features of SNe la, although it cannot be taken as a 
confirmation of the hypothesis, leaves the possibility of a flame- 
driven DDT open to a more detailed investigation of potential 
physical mechanisms. This result is not obvious, since the time 
and location of the DDT is not fixed and it thus may occur at 
instances disfavoring an efficient detonation. However, the three 
simulations presented here trigger a successful detonation after 
very different deflagration phases indicating robustness of the 
criterion. Parameterizing the DDT according to properties of the 
turbulent deflagration flame is naturally very sensitive to the tur- 
bulence description and we emphasize the importance of a well- 
founded subgrid-scale model. 

The three-dimensional simulations presented here produce 
somewhat less energetic explosions than the two-dimensional 
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Table 1. Model parameters: the DDT columns provide values for the time of the initiation of the detonation wave (JddtX f° r its 
distance from the center of the WD, for the density at the DDT spot (pddtX an d f° r the nuclear energy release in the deflagration 
prior to the DDT (.E™ 1 ). The last six columns characterize the asymptotic kinetic energy of the ejecta (E^ mpt ) and the composition 
of the final stages of the explosion models, and give the ratio Wjff of the masses of nuclear ash produced by the deflagration to that 
produced by the detonation. 
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center [10 s cm] 


[10 7 g 


[foe] 


[foe] 


[M ] 


[M a ] 


[M Q ] 


[M ] 


DD_005 


5 


0.731 


1.733 




1.33 


0.221 


1.524 


1.141 


0.220 


0.00586 


0.0386 


0.37 


DD.020 


20 


0.724 


1.733 




1.92 


0.500 


1.237 


0.833 


0.435 


0.0190 


0.118 


1.48 


DD_800 


800 


0.675 


1.820 




2.40 


0.656 


1.004 


0.638 


0.547 


0.0297 


0.189 


1.59 



setups of Golombek & Niemever (2005j). This is most likely due 
to the symmetry constraints in the two-dimensional simulations 
where the detonation wave is initiated in two toroidal rings. 
Moreover, deflagrations are more energetic in three dimensions 
and produce a more complex structure do to the additional de- 
gree of freedom, both effects hampering the detonation wave 
propagation. Nonetheless, almost all WD material is subject to 
nuclear burning in the models studied here. Yet the detonation 
wave does not in all cases outrun the deflagration features on the 
far side of the DDT. This leads to asymmetries in the structure 
of the outermost ejecta which, however, contain little mass. Such 
asymmetry effects may be an artifact of our restriction to a single 
DDT. The DDT criterion can easily be modified to account for 
multiple DDTs at each patch of the flame where the distributed 
burning regime is entered. Such parameters of the DDT criterion 
will certainly have an impact on the explosion which will be ana- 
lyzed in a forthcoming study. Although we expect that the range 
of variation in the results will be similar to that found with a sin- 
gle DDT, the correlation between the number of ignition kernels 
and the explosion strength will possibly shift towards stronger 
explosions for denser ignitions. 

The results of the explosion simulations are in good agree- 
ment with current observational constraints of SNe la. Our set 
of models indicates that the presented delayed detonation sce- 
nario can cover at least large parts of the observed sample of 
SNe la with respect to explosion energies and NSE material pro- 
duction. However, for a thorough comparison with observational 
results, a detailed analysis of the ejecta composition on the basis 
of high-resolution models is required and will be addressed in a 
forthcoming study. 
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